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Abstract 

We study the perturbative QCD corrections to heavy-quark structure functions of charged-lepton 
deep-inelastic scattering and their impact on global fits of parton distributions. We include the log- 
arithmically enhanced terms near threshold due to soft gluon resummation in the QCD corrections 
at next-to-next-to-leading order. We demonstrate that this approximation is sufficient to describe 
the available HERA data in most parts of the kinematic region. The threshold-enhanced next-to- 
next-to-leading order corrections improve the agreement between predictions based on global fits 
of the parton distribution functions and the HERA collider data even in the small- jc region. 
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In this letter we study structure functions in deep-inelastic scattering (DIS) of charged leptons off 
nucleons. We focus on the production of heavy quarks, like e.g., charm, which proceeds within 
perturbative Quantum Chromodynamics (QCD) predominantly through boson-gluon fusion [1,2]. 
This is a reaction of great interest because at moderate momentum transfer Q, it provides a direct 
probe of the gluon content of the nucleon over a wide range of Bjorken x. In the present work, 
we specifically include higher order QCD corrections to DIS heavy-quark production and we wish 
to determine their impact on our information about the parton distribution functions (PDFs) and, 
especially, the gluon PDF. 

Our motivation stems from the high statistics data for the charm structure function provided 
by the HERA experiments [3,4], where accounts for a large fraction (up to 30%) of the total 
structure function F2. The presently available DIS data allows for high precision extractions of 
PDFs in global fits [5,6]. The treatment of the charm contribution in these fits is an important issue 
as it can induce potentially large effects also in the PDFs of light quarks and the gluon obtained 
from these global fits (see e.g. the recent review [7]). 

In the standard factorization approach the heavy-quark contribution to the DIS structure func- 
tions can be written as a convolution of PDFs and coefficient functions, 
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Fk{x,Q\m^) = ^^^ ^ JdzA(z,iuj)ci,k(r]{x/z),^,luj,iu^), (1) 



where a = I +4m /Q and m, eq are the heavy-quark mass and charge. The strong coupling 
constant at the renormalization scale yu^ is denoted as = (Xsip-r) the PDFs for the parton of flavor 
/ are fi(x,/up at the factorization scale yu^. Both scales //^ and /i,. are assumed equal throughout 
this work, i.e. /i = //^ = /i^. The kinematical variables rj and ^ in Eq. ([T]) are given as 

_ s ^ 

and the partonic center-of-mass energy s = Q^iz/x- 1) in Eq. ([T]). The coefficient functions of the 
hard partonic scattering process enjoy an expansion in as as 

00 k 2 

Ci,k(r],^,iu^) = _2(47ra,0u))^ ^cjf (77,^)ln'^. (3) 
k=o 1=0 ' 

The perturbative QCD predictions for the coefficient functions to the leading order (LO) are 
long known [1,2]. The next-to-leading order (NLO) radiative corrections are available since more 
than 15 years [8]. Likewise, the massless coefficient functions for the light-quark content of the 
DIS structure functions have also been calculated to the next-to-next-to-leading order (NNLO) 
some time ago [9-13]. More recently, the scale evolution of the PDFs has been matching the 
NNLO accuracy [14,15]. 

Although the full heavy-quark coefficient functions at two loops are the only unknown for a 
complete NNLO analysis, much can already be said about these terms at present. Because the 
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structure functions for massive quarks contain two hard scales, that is the momentum transfer Q 
and the heavy-quark mass m, the study of particular kinematical limits yields valuable information. 

For instance, at asymptotic values Q^,n^ » one may treat the heavy quark as effectively 
massless. As an upshot, large logarithms Xnijj.^ /np-) are summed over in the PDF evolution [16]. 
Accordingly, the higher order coefficient functions [8] in Eq. ([3]) take asymptotic forms [17, 18] 
and, together with the PDFs, require matching. This is the standard procedure when changing 
the description from QCD with light flavors and a heavy quark to a theory with n^r + l light 
quarks. Thus, a so-called variable flavor number scheme (VFNS) has to describe this transition in 
the effective number of light flavors and, moreover, a general-mass formalism for a VFNS has to 
be consistent with QCD factorization, see Refs. [19,20]. However, the VFNS formalism cannot 
be routinely extrapolated to the region of 2 ~ m, where the power corrections due to the heavy- 
quark mass effects break the factorization and therefore the massive quarks cannot be considered 
as partonic constituents of the nucleon. In this kinematical region, one strictly applies QCD with 
Uj: light flavors and one heavy quark, thus working in the so-called fixed flavor number scheme 
(FFNS). The VFNS can be used for the whole kinematics of the existing DIS data only if it is 
matched in order to provide a smooth transition to the FFNS at small Q. Details of this transition 
cannot be derived within the VFNS framework and are therefore subject to model assumptions 
(see Ref. [7] for the review of the modem state of art and history of this modeling). On the other 
hand, in the FFNS the large logarithms Xnijj.^ I irp-) are contained in the higher order corrections to 
the coefficient functions. Thus, the basic motivation for the use of a VFNS weakens, at least for 
the realistic DIS kinematics [21]. 

Near threshold, for s - 4m^ or equivalently ?7 «; 1, higher order perturbative corrections are 
much enhanced. The coefficient functions in Eq. ([3]) exhibit large double logarithms or /ln^'yS, with 
jS = sjl -AnP-/ s being the velocity of the heavy quark and these Sudakov logarithms can be re- 
summed to all orders in perturbation theory. Currently this has been achieved to the next-to-leading 
logarithmic (NLL) accuracy and the resummed result can be employed to generate approximate 
results at NNLO in QCD (see e.g. [22]). In the present work, we focus on the impact of these 
approximate NNLO corrections due to soft gluons and assess their impact on global fits of PDFs 
in the FFNS. To that end, it is instructive to express the convolution in Eq. ([I]) as an integration 
over the partonic variable 77. In this way, we obtain with z{rf) = .x:(l -1- 4(7/ -i- 1)/^) 



where the integration is bounded by rimax = ^/4(1 /x-l)-l. 

In Fig. [H we plot the shape of the gluon distribution g(z) = fg{z) for representative kinematics of 
charm quark production at the HERA collider. The typical values of x and Q employed for xg[zirj)] 
in Fig. [T] are correlated: The minimal (maximal) value of x corresponds to the minimal (maximal) 
value of Q. For small scales Q the value of xg[z{ri)] is suppressed at large rj due to its argument, 
because z(rj) rises with rj about linearly (see Eq. dH)) and the gluon PDF decreases with rising 
argument. Therefore, for small scales Q the region of 77 ^ 1 provides the dominant contribution 
to the charm structure function and the parton kinematics relevant for the coefficient functions 
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Figure 1 : The Ty-dependence of the gluon distribution xg{z) for representative kinematics of the HERA 
collider experiments (see Eq. dUl). 



are effectively constrained to the threshold region (see Ref. [23]). Recall that is dominated by 
boson-gluon fusion. On the other hand at large virtualities Q the rise of z(ri) with r] is not so fast. As 
a consequence, the suppression of the large-// region due to the shape of the gluon PDF is weaker, 
as can be seen in Fig.[Il However, for the case of bottom-quark production the large-?/ suppression 
is stronger due to the larger quark mass. Hence, even at large scales Q the bottom structure function 
is still saturated by parton processes (and coefficient functions) close to threshold. 

Let us now briefly summarize the threshold approximation to the hard scattering coefficient 
functions in Eq. (|3]). To that end, we are following the standard procedure for resumming Sudakov 
logarithms, see e.g. Refs. [24,25]. In differential kinematics for the one-particle inclusive DIS 
production of a heavy quark (see [8,22]) the (dominant) gluon coefficient function is given by 

/(l+/?)/2 -C" ,2 (0,0), , , , 

am C C r^ d^c\ \s ,ti,ui) 

= J ^(-^i) J ds4 K^'\s',h,ui) , (5) 

i'(l-/^)/2 

where cl'^''^^ is the Bom contribution. We have ui = s' + t\- and s' = s + Q^ and 

In a physical interpretation denotes the additional energy carried away by soft gluon emission 
above the partonic threshold. At higher orders, Eq. ^ contains plus-distributions of the type 
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a[[\n^'~^(s4/m^)/ S4]+ that give rise to the Sudakov logarithms upon integration, i.e. the well 
known double logarithms aJ\n^'/3 of an inclusive formulation. At the differential level (one- 
particle inclusive kinematics) the threshold resummation for DIS heavy-quark production has been 
performed to NLL accuracy in Ref. [22]. Subsequently, the resummed result has been used to 
generate the factors K^'^ (cf. Eq. ([5])) at fixed-order perturbation theory through NNLO. These 
factors K^'^ contain the large logarithms. 

In the present paper we improve the approximate NNLO results of Ref. [22] by performing a 
matching at one-loop and by including the NLO Coulomb corrections. This provides us the with 
the first three powers of Sudakov logarithms at all orders and we arrive at the following expressions 
through NNLO, 
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Cj,{Lp + \n{hlui))-lCp{\+Lp) 
+21n(l - r,2)L^ + ln(r,)L^) + \cj^n{ti Imf- 3^2 - 41n(l - r,^)Lp 
+2 ln(r,)Ly3 - 2 ln(r,) \n{ti /w i ) - ln(r,)^ + 2Li2( 1 -ui/ti/rs) + 2Li2( 1 -h/ui /r,) j 
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C^K + 2C^Cp(l- 3Lj3 - 2L2 - 2 Inih /w 1 ) - 2 Inih /m 1 )Ly3 + 2 ln( 1 - r/)Ly3 

+ ln(r,)Lpj + CpM^ + Lp) + ACj{l + Lpf - ^Cj^pJ^Lp + In^i /wo) 

+ ^ Cj(2L^ + 4 InC?! /m 1 )L/3 + 3 InC?! /m 1 f + 4 ln(-M 1 /m^) - 1 1 ^2 - 4 ln( 1 - r,2)L/j 
+2\n{rs)Lp - 21n(r,)ln(fi/Mi) -ln(r,f + 2Li2(l -Mi/?i/r,) + 2Li2(l - ^i/wi/r,)) 
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(^2-2Li2(-r,)-2Li2(r,)) 



where D/ = [In (54/m )/54]+ denote the plus-distribution. We have in QCD = 3, = 4/3, 
jSo = 11/3C^ - 2/3ny and K = (67/18 -^2)^^ - 5/9n^. The variables and and Lp are given by 
r, = (1 -/3)Kl +/3) and Lp = (1 -2m^/s)//3 {ln(r,) + i7r}. 

Our improved NNLO approximation in Eq. ^ is exact in the region of phase space s - 4m^, 
where perturbative corrections receive the largest weight from the convolution with the gluon PDF 
(see discussion above and Refs. [22, 23]). In Eq. ^ we have restricted ourselves to the case 
yu^ = m^. While it is straightforward to allow for general choices //^ ^ to logarithmic accuracy 
in the threshold resummation formalism one can even derive the exact /u^ and //^^ scale dependence 

through NNLO [22] with the help of renormalization group methods. Thus, the functions C2 ^ and 



c?'^^ in Eq. ^ are known exactly [22], see Fig.|2]for plots. 
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Figure 2: (a) The ?7-dependence of the gluon coefficient functions at NNLO contributing to the heavy- 
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quark DIS structure function F2 (cf. Eq. dS])). The solid Hne denotes c? '^^ according to Eq. ^ dashes and 



(2 1) (2 2) i1 0) 

dots show the exact result for c\ , and ' from [22]. (b) The low-x (large-77) asymptotics of cI, of 

J V J 4'^ J L 

Eq. (1101) (solid line) from [26] and Eq. (1121) as modeled in [27] (dashes). Eq. (1121 ) without account of the 
model suppression factor (dashed dots) and the model subtraction term (dots) are given for illustration. 



In a dilferent kinematical regime, small-x effects in DIS heavy-flavor production have been 
studied systematically [26] and incorporated in phenomenological analyses [27]. Based on kj- 
resummation [26], the leading logarithm at small-x for c?'*^^ can be derived: 

4?('^'^) = Aw-)^' (10) 

recall z(ri) = x(l+ 4(?7 +1)/^). The function K2 is a low order polynomial in ^ = Q^/rrP- that can be 
determined from empirical fits to the leading ln(l /x) behavior of Ref. [26]. In the region of ^ ^ 1 
it has been estimated as [27] 

y^2(^) = 13.073^-23.827^2 + 24.107^3_9 i73^_ (11) 

However, for phenomenological applications the sole knowledge of ln(.x;)-terms is usually insuf- 
ficient and additional assumption have to be supplied. Thus, a particular functional form for the 
coefficient function at small-x has been suggested in Ref. [27], 

^f^'?'^) = TT^f^ (Mz/x)-A) (1 -ax/z)20 ^ (12) 

where the subtraction term (-4) at ln(l/z) is motivated by the terms of similar size in the small- 
X limit of other known coefficient and splitting functions. The factor yS times the polynomial 
(1 - ax/zf'^ suppresses large-z effects by a large power and a is given below Eq. ([T]). 
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We plot the small-x asymptotics of ci in Fig. [21 where we display Eq. (fTOl) for the original 
result of Ref. [26] and the model coefficient function of Eq. (1121) . The variants of Eq. (1121) without 
account of the subtraction term (-4) and the suppression factor of /3(\ -ax/zf''^ are also given for 
comparison. The form of Eq. (fT2l) is very sensitive to the particular choice of these terms. Due 
to the large-z suppression Eq. (fT2l) vanishes at rj < 20, and the region of 77, where this happens 
is defined by the power of 20 of the polynomial. Similarly, at large 77 the form of Eq. (fT2l) is 
entirely defined by the subtraction term. The presence of subleading terms in the small-x expansion 
with large numerical coefficients is a well documented feature at higher orders in QCD, see e.g. 
the example of the three-loop gluon splitting function [15]. Moreover, precise phenomenological 
predictions based on small-x approximations only are extremely difficult to make, because the 
convolution in Eq. (U) is non-local. Thus, the small-.x; terms in the coefficient function are weighted 
by the PDFs in the large- jc region (and vice versa), see [15]. 
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Figure 3: The predictions of the global fit of PDFs for the charm structure function F2 compared to data 
of Ref. [4] on obtained in Run I of the HERA collider. The QCD corrections to the charm coefficient 
function in Eq. ([3]) have been included up to NNLO (solid lines), NLO (dashes), and LO (dashed-dots) and 
the scale has been chosen fi^ = +4nP-. The variant of the NNLO fit with the scale choice /i = 2 is given 
by dots. 

For phenomenological applications our approximate NNLO result Eq. ^ for c?'^'' is added 
on top of the exact NLO predictions [8] and supplemented with the exact NNLO scale dependent 
functions c^2g'' '^2^g^ of Ref. [22]. This provides the best present estimate for, say, the nucleon 
structure function for DIS electro-production of charm quarks. We investigate the impact of 
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these NNLO (gluon induced) contributions to on the nucleon PDFs extracted from global fits 
and perform a modified version of the fit of Ref. [6]. That fit is based on global data on inclusive 
charged-lepton DIS off" nucleons supplemented by data for dimuon nucleon-nucleon production 
(i.e. the Drell-Yan process). Within the FFNS we take into account the NNLO corrections to the 
QCD evolution [14, 15] and the massless DIS and Drell-Yan coefficient functions [9-12,28-30]. 
The value of m was fixed at 1.25GeV, close to the world average and the factorization scale was 
selected as sJQ-^ + 4m-^. For comparision we also provide NNLO results for the factorization scale 
a.tfi = Q. The latter choice naturally leads to larger deviations at small virtualities Q due to a much 
increased numerical value of as (note that we have identified ju= iUj:= ju^). 

Q^=2 GeV^ 
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Figure 4: The impact of the large-?? tail of the NNLO coefRcient functions on the value of at small Q 
displaying our calculation (solid curve) and the same with the model of Ref. [27] for Cj^^"^ at large-77 added 
(dashes). The NLO calculation (dashed-dots) and the data of Ref. [4] are given for comparison. 




We also perform two variants of this fit taking into account only the LO corrections of Refs. [1, 
2] and the NLO corrections of Ref. [8]. The predictions for F2 based on these three fits are 
compared in Fig. [3] to the ZEUS data of Ref. [4]. The latter data are not used in the fits. At 
the smallest values of x and Q in the plot the predictions rise monotonically with increasing orders 
of perturbative QCD, thus improving agreement with the data. As we discussed above, in this 
region the value of F2 is not sensitive to the coefficient functions at large 77 and therefore our 
predictions (cf. Eq. ^) can be considered as a good approximation to the full NNLO result for 

(2 I) 

F^- At bigger values of x and Q, as a result of a negative contribution from c^^ at large 77 the 

NNLO predictions dip below the NLO ones. In this region of rj the value of c?'*^^ was set to 
zero as our choice of matching the threshold approximation to fixed order perturbation theory (see 
e.g. [31] for related discussions). Checking the curves of Fig. [3] at large values of x and Q one can 
conclude that this contribution should be positive in order to improve the agreement with the data, 
the particular numerical impact depending, of course, on the gluon distribution shape as one can 
conclude from Fig.[T] In this region (large x and Q), the slope of F^ appears to be distinctly ffatter 
in Q, particularly at higher x. However, the kinematics for large values of Q is far from threshold 
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and beyond control of our soft gluon approximation. 

The small-x contribution to c^2g'' modeled in Ref. [27] affects the comparison to data at the 
lowest Q and at small-;c only. In order to assess the impact of the small-x term of Eq. (flOl) quan- 
titatively, we focus on the lowest bin = 2 GeV and illustrate its effect in Fig.lH At the lowest 
value in x we do observe a slight sensitivity on the small-x term, which in terms of t]{z) corresponds 
to the region of larger 77 » 1 . The effect of the high-;7 model of Ref. [27] amounts at most to a 
30%-fraction of the dominant NNLO contribution coming from the threshold region at the lowest 
values of x. However, in this region, the addition of the small-x contribution to c?'*^^ overshoots 

6 ' 2,g 

the data while it vanishes quickly at larger x. Recall in our analysis we set the contribution of the 
threshold logarithms in c?'*^'' to zero for 77 > 1 . As we discussed above, the ansatz of Eq. (flOl) has 
an inherent model uncertainty of 100% since it is driven by ad hoc parameters. Therefore it cannot 
be used in quantitative comparisons. 
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Figure 5: The same as Fig. |3]for the charm electro-production cross section data obtained by the EMC 
experiment [32]. 

At fixed-target energies the charm contribution to the inclusive sample is much smaller than 
at the HERA collider. This makes the experimental determination of F2 more difficult. The only 
conclusive data on fixed-target charm electro-production were obtained some time ago by the EMC 
collaboration [32]. These data are compared to the predictions of our fits in Fig. [51 At smallest 
values of Q the data are in disagreement with the predictions. Moreover the disagreement increases 
with the order of the perturbative QCD correction. At bigger values of Q the difference between 
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the NLO and the NNLO predictions is marginal due to the difference in the gluon PDFs obtained 
in these variants of the fit. However the general agreement between data and calculations is far 
from ideal. Due to the limited collision energy the EMC data are sensitive to the region of ?/ < 1 
only, where our NNLO approximation in Eq. Q should describe the exact coefficient function 
^^ig"* well. Thus, we see no way to improve the agreement with the EMC data by performing 

a complete calculation of c^2g^- Since the EMC data are unique it seems to be useful to have 
additional experimental input to clarify this disagreement. 
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Figure 6: The ±\cr band for the gluon (left panel) and non-strange sea (right panel) distributions obtained 
in the fit by including the coefficient functions for charm electro-production up to NNLO (solid lines) and 
NLO (dashes). 



Despite the fact that we do not use data for the charm electro-production in the fit our results 
are sensitive to the details of the description of F^. This is because F'j amounts to a substantial 
contribution, up to 30%, to the inclusive DIS structure functions at the HERA collider energies. 
The biggest variation in the fitted PDFs due to the NNLO corrections to is observed for the sea 
quarks PDF at small-x. The latter becomes smaller in order to compensate the positive contribution 
from the NNLO term at small-x (see Fig. (6]). The gluon distribution at x ~ 0.02 becomes bigger 
by about one standard deviation as a compensation of the negative contribution of the NNLO term 
and the fitted value of as(Mz), which is anti-correlated with the gluon distribution at small-x, goes 
down by about Icr. Other PDFs are essentially not affected by the corrections. 

Let us summarize: We have improved the perturbative QCD predictions for the heavy-quark 
DIS structure functions. Our NNLO approximation takes along the first three powers of Sudakov 
logarithms for the boson-gluon channel yg, which, in an inclusive formulation (performing the 
integration and keeping the leading terms in /3 only) corresponds to all logarithmically enhanced 
terms In^yS, k = 2,...,4-. Moreover, we have employed the exact expressions for all scale dependent 
terms through NNLO [22]. 
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Subsequently, we have applied the NLO QCD corrections [8] to the charm structure function 
and our new approximate NNLO result in a global fit to data for charged-lepton DIS and dimuon 
production in the Drell-Yan process. Especially the use of the threshold-approximated NNLO re- 
sult is legitimate because the gluon PDF constrains the parton kinematics to values around s - 4m^. 
Modifying the fit of Ref. [6] in this way we have studied the eff"ects for the determination of PDFs. 
We have found that our approximate NNLO result for gives better agreement between the fitted 
PDFs and the HERA collider data. The agreement with data extends even down to small values 
of X. The results of our fit are also in good agreement with ZEUS data [4] on charm electro- 
production, which was not used in the fits. Comparing with EMC data [32] we did find disagree- 
ment, though, and it would be interesting to get new and independent experimental information in 
order to resolve it. 

On the theory side we could in principle extend the NNLO threshold approximation of Eq. dH) 
further to include (after integration) the linear logarithm in ln/3 at two loops and the two-loop 
Coulomb corrections following the procedure of Ref. [31] for heavy-quark hadro-production. 
However, we leave this to future research. 
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